Faithful Effective-One-Body waveforms of 
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We compare a recently derived, resummed high post-Newtonian accuracy Effective-One-Body (EOB) 
quadrupolar waveform to the results of a numerical simulation of the inspired and merger of an equal-mass 
black hole binary. We find a remarkable agreement, both in phase and in amplitude, with a maximal dephasing 
which can be reduced below ±0.005 gravitational-wave (GW) cycles over 12 GW cycles corresponding to the 
end of the inspiral, the plunge, the merger and the beginning of the ringdown. This level of agreement is shown 
for two different values of the effective fourth post-Newtonian parameter a 5 , and for corresponding, appropri- 
ately "flexed" values of the radiation-reaction resummation parameter v po \ e . In addition, our resummed EOB 
amplitude agrees to better than the 1% level with the numerical-relativity one up to the late inspiral. These 
results, together with other recent work on the EOB-numerical-relativity comparison, confirm the ability of the 
EOB formalism to accurately capture the general relativistic waveforms. 

PACS numbers: 04.25.Nx, 04.30.-w, 04.30.Db 
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I. INTRODUCTION 



The gravitational-wave (GW) signals emitted by coalesc- 
ing black hole binaries are among the most promising targets 
for the currently operating network of ground-based detectors 
GEO/LIGO/VIRGO. The most useful part of the waveform 
for detection comes from the most relativistic part of the dy- 
namics, around the coalescence, i.e. the last few cycles of 
the adiabatic inspiral, the plunge and the merger. It is crucial 
for GW detection purposes to have available a large bank of 
"templates" that accurately represent the GW signals radiated 
by these binaries. The construction of faithful 1 GW templates 
for coalescing binaries comprising spinning black holes (with 
arbitrary masses mi, m.2 and spins Si, S2) is a non trivial 
task. In view of the multi-dimensionality of the corresponding 
parameter space, state-of-the-art numerical simulations can- 
not densely sample this parameter space. It is therefore ur- 
gent to devise analytical methods for computing (as a function 
of the physical parameters mi, r«2, Si, S2) the correspond- 
ing GW waveforms. Here we continue the program of con- 
structing, within the Effective-One-Body (EOB) method [2- 
5] high-accuracy analytic waveforms describing the GW sig- 
nal emitted by inspiralling and merging binary black holes 
with arbitrary masses and spins. The EOB method was the 
first to provide estimates of the complete waveform (cover- 
ing inspiral, plunge, merger and ring-down) of a coalescing 
black hole binary, both for non-spinning systems [3], and for 



spinning ones [6]. 

Numerical Relativity (NR) recently succeeded in giving us 
access to reliable information about the dynamics and radia- 
tion of binary black hole coalescences [7-18]. This opens the 
possibility of comparing the EOB predictions to NR results. 

The comparison between the EOB approach and NR results 
has been recently initiated in several works [19-24]. These re- 
cent comparisons have been done using two different versions 
of EOB waveforms. The works of Buonannoet al. [19, 20, 23] 
used a restricted waveform, as proposed in the first EOB pa- 
per [3], but with an improved matching to the ringdown (sim- 
ilar to the one used in [25]) making use of three quasi-normal 
modes. By contrast, the recent works of Damour and Na- 
gar [22, 24] use a new, resummed high post-Newtonian (PN) 
accuracy 2 EOB quadrupolar waveform. This improved EOB 
waveform has been shown to exhibit a remarkable agreement, 
both in phase and in amplitude, with NR waveforms in two 
separate physical situations: (i) inspiral and coalescence of 
small-mass-ratio (non-spinning) systems [22] (comparing it 
to waveforms computed by means of numerical simulations 
of test particles, with an added radiation-reaction force, mov- 
ing in black-hole backgrounds [26]) and (ii) inspiral (up to a 
limiting GW frequency ~ 0.14/M) of an equal-mass (non- 
spinning) system [24] (comparing it to recently published re- 
sults of a high-accuracy inspiral simulation [18]). 

The present paper is a continuation of the general program 
of constructing, within the Effective-One-Body (EOB) ap- 
proach, high accuracy, faithful analytic waveforms describing 



Following the terminology of [1], we recall that "effectual templates" are 
templates exhibiting large overlaps with an exact signal after maximizing 
over all (kinematical and dynamical) parameters, while "faithful" ones are 
so "close" to an exact one that they have large overlaps for values of the dy- 
namical parameters which are very close to the real ones ("small biases"). 



2 This high PN accuracy can be called 3 +2 -PN because it includes not only 
the known comparable-mass 3 PN waveform corrections, but also the test- 
mass limit of the 4 PN and 5 PN waveform amplitude corrections. See [24] 
for details and references. 
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the gravitational wave signal emitted by inspiralling and co- 
alescing binary black holes. Here we shall consider the coa- 
lescence signal emitted by a non-spinning equal-mass binary 
black-hole system. We shall compare the phase and the ampli- 
tude of the new resummed 3 +2 PN-accurate EOB quadrupolar 
waveform of [22, 24] to a numerical relativity simulation of a 
coalescing black hole binary performed at the Albert Einstein 
Institute (AEI). 

This comparison will confirm the ability of the EOB ap- 
proach to provide accurate analytical representations of NR 
waveforms. We note that the recent work [23] had already 
shown the ability of the analytically less accurate restricted 3 
EOB waveforms to provide rather accurate approximations to 
NASA-Goddard NR coalescence waveforms for several dif- 
ferent mass ratios (mi/m2 = 1,3/2,2 and 4). More pre- 
cisely, Ref. [23] found, in the equal mass case, an EOB/NR 
dephasing of ~ ±0.03 GW cycles over 15 GW cycles. Con- 
cerning the amplitude, the latter reference does not quan- 
tify the restricted EOB/NR difference, but one can read from 
Fig. 21 of [18] that the difference between the restricted PN 
(or EOB) amplitude and the Caltech-Cornell inspiral NR one 
is ~ 7%. By contrast, the present paper will show that the 
new, resummed waveform exhibits a significantly smaller de- 
phasing ~ ±0.005 GW cycles over 12 GW cycles, and, most 
remarkably, exhibits an excellent agreement in amplitude, 
both during the inspiral (where it is better than the 1% level) 
and the ring-down. This good result is obtained by making 
use (as proposed in several previous works [5, 23, 24, 27, 28]) 
of the natural flexibility of the EOB approach. 

An alternative approach to the construction of analytical 
templates to model (non-spinning) coalescing binary black 
holes with arbitrary mass ratios has been recently proposed 
in Refs. [29,30]. 

This paper is organized as follows: In Sec. II we briefly 
describe the numerical simulation, whose results we use in 
the following. In Sec. Ill we spell out the features of the EOB 
waveform that we shall use. The main section is Sec. IV where 
we compare the new, resummed EOB waveform to NR data. 
We also include a comparison where we use the less accurate 
"restricted" EOB waveform, and simpler QNM-matching, 
used in some of the previous EOB works [3, 19, 20, 23]. The 
paper ends with some conclusions. 



n. BRIEF DESCRIPTION OF THE NUMERICAL 
SIMULATION 

The numerical simulations have been carried out with the 
Ccatie code [31], a three-dimensional finite-differencing 
code developed at the Albert Einstein Institute and at the Cen- 
ter for Computation and Technology (CCT) of the Louisiana 
State University. The code is based on the Cactus Compu- 
tational Toolkit [32] for the solution of the Einstein equations 
in a finite-size domain covered with a Cartesian rectangular 



grid. The main and new features of the code have been re- 
cently discussed in Ref. [31], and we here briefly recall the 
most important ones only. 

The Einstein equations are formulated as an initial-value 
problem via a conformal and traceless "3 + 1" decomposi- 
tion. The spacetime geometry is decomposed into: (i) the 3- 
metric of spacelike slices, (ii) the extrinsic curvature of those 
slices, and (iii) the lapse and shift. See [31] for the explicit 
form of the equations. The lapse function is evolved using the 
"1 + log" slicing condition [33], while the shift is evolved us- 
ing the hyperbolic T-driver condition discussed in Ref. [34], 
but with the difference that advection terms have been added 
following the experience of [8, 35], and are required for cor- 
rect advection of the punctures in "moving-puncture" evolu- 
tions. 

Spatial differentiation of the evolution variables is per- 
formed via straightforward finite-differencing using fourth- 
order accurate centered stencils for all but the advection terms 
for each variable, which are instead upwinded in the direction 
of the shift. Vertex-centered adaptive mesh-refinement is em- 
ployed using nested grids via the Carpet infrastructure [36], 
with a 2 : 1 refinement for successive grid levels, and the 
highest resolution concentrated in the neighborhood of the in- 
dividual horizons. Individual apparent horizons are located 
every few time steps during the time evolution [37], which is 
obtained via a "method-of-lines" and with a fourth-order ac- 
curate Runge-Kutta time integrator. 

The simulations were performed on a domain with outer 
boundaries located at 4 768M C , and a grid structure consists 
of nine mesh-refinement levels, the finest of which has a spa- 
tial resolution of h = 0.02M C . Simulations with lower res- 
olution (i.e., with h = 0.024M C and h = 0.03M C ) have 
also been carried out to validate the consistency of the re- 
sults. An important feature of the Ccatie code is the pos- 
sibility of employing two distinct methods for the calculation 
of the gravitational radiation produced. The first method uses 
the Newman-Penrose curvature scalar ip^, with respect to a 
suitable frame at the extraction radius. An alternative method 
measures the metric of the numerically generated spacetime 
against a fixed background at the extraction radius, and de- 
termines the gauge-invariant Regge-Wheeler-Zerilli-Moncrief 
functions (see Ref. [38] for a review and references). Both 
methods have been systematically studied in Ref. [31], where 
they were also compared and shown to yield essentially iden- 
tical results, both in terms of their asymptotic scaling prop- 
erties (e.g., the peeling-theorem), and in terms of the polar- 
ization amplitudes h+ and h x . The analysis carried out here 
used as basic NR data the gauge-invariant (Zerilli-Moncrief) 
metric perturbations. These were extracted on (NR) coordi- 
nate 2-spheres with (NR) coordinate radii i?NR = 60M C up 
to i?NR = 120M C , with a separation of 10M C between two ad- 
jacent observers. The analysis carried out below uses, as ap- 



3 Here "restricted" refers to a waveform h oc Q 2 / 3 e 2l *. 



4 We denote by M c the internal length and mass units used in the code (with 
G = c = 1). Beware that M c slightly differs from M = m\ + ni2 (see 
below). 
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proximate asymptotic amplitude, the metric perturbation ex- 
tracted at i? NR = 120M C . 

The initial data for the black-hole binary are obtained 
by a Brill-Lindquist [39] construction, where the additional 
asymptotically flat end of each wormhole is compactified into 
a single point, the so called puncture [40]. This approach ex- 
plicitly uses the Bowen-York extrinsic curvature and solves 
the Hamiltonian constraint equation numerically (as detailed 
in Ref. [41]), after having chosen the free parameters for 
the puncture initial data. Quasi-circularity of the initial orbit 
can then be obtained by specifying the puncture parameters 
in terms of an effective-potential method [42] as discussed 
in [31]. However, the assumption of "quasi-circularity" (in 
the sense of [42]) at the (rather small) initial separations fre- 
quently used in numerical-relativity simulations has the draw- 
back of introducing a small but nonzero amount of eccentric- 
ity. To compensate for, or reduce, this effect, other approaches 
have been suggested recently. One of these is based on an 
iterative minimization procedure where, throughout a series 
of simulations with slightly different initial black hole con- 
figurations, the eccentricity is measured and minimized [43]. 
A simpler and rather effective approach has been proposed 
in Ref. [44], and consists of specifying the initial puncture- 
parameters as the end-state of a binary system whose evolu- 
tion is determined, starting from a large separation, via the 
solution of the Taylor-expanded 3 PN-accurate equations of 
motion [6, 45, 46]. 

We have here essentially followed this latter prescription 
and considered, in particular, the initial data denoted by Ell 
in Table I of [44], that have been shown there to reduce the 
eccentricity to e < 0.002. More specifically, our initial black 
holes have a coordinate distance D = 11M C , momenta in the 
radial and tangential directions of P r = -7.09412 x 10~ 4 Af c 
and P t = 0.0900993M C , and a puncture mass-parameter of 
0.487035M C , leading to initial individual black-hole masses 
mi = m,2 = 0.499821M C , and thus a total mass of the binary 
system M = mi + mi = 0.999642M C . Overall, the sim- 
ulation covers about ~ 1600 M of the final evolution of the 
binary, thus comprising 8 orbits and about 16 GW cycles. 

The mass and spin of the final black hole have been com- 
puted through two different methods yielding, however, very 
similar results: (i) by using the isolated/dynamical horizon 
formalism [47, 48], where a proper rotational Killing vector is 
searched on the final apparent horizon to measure the spin, 
and the horizon area is used for computing the black hole 
mass (see Sec. IV D of Ref. [31] for details); (ii) by perform- 
ing a fit of the dominant quasi-normal mode 5 of the complex 
ringdown waveform. This fit was performed by a non-linear 
least-squares Gauss-Newton method, using exp(— at + p) as a 
parameter-dependent template (with two complex parameters 
(a, p)), and an appropriate time interval during the ringdown 
(chosen by minimizing the post-fit residual). [For a discus- 
sion of methods for QNM fitting see Refs. [49-51]]. Then, 



5 In the notation introduced in Sec. Ill below, the dominant mode corre- 
sponds to the labels (±, £, £', m, n) = (+, 2, 2, 2, 0). 



from the best-fit value of a (i.e., the QNM dominant complex 
frequency 03220)' we computed the values of the mass and di- 
mensionless spin parameters of the final black hole by using 
the interpolating fits given in Appendix E of Ref. [52]. The 
results of these two methods are denoted as (M hor , j hor ) and 
(M rin s, j rin s), respectively. 

The most relevant properties of the binary system are sum- 
marized in Table I. The difference (which is < 1%) be- 
tween the quoted values of the final black hole parameters 
might come, in part, from inaccuracies in the interpolating 
fits of Ref. [52]. In the following we will use, in our EOB- 
matching procedure, the ringdown-fitted black hole parame- 
ters (M rmg , j rmg ) (so that the dominant complex frequency 
will be guaranteed to have the best possible value). 



III. EFFECTIVE-ONE-BODY (EOB) METHOD AND 
WAVEFORM 

We shall not review here in detail the EOB method [2-5], 
which has been described in several recent publications, no- 
tably Refs. [23, 24]. We shall only indicate the EOB elements 
that are crucial for the present study. For detailed definitions 
of the EOB ingredients we refer to the recent paper [24] that 
we follow, except when otherwise indicated below. 

Before entering the details of our EOB implementation, let 
us recall that Ref. [24] proposed a methodology for improving 
the waveform implementation of the EOB philosophy based 
on understanding, element by element, the physics behind 
each feature of the waveform, and on systematically compar- 
ing various EOB-based waveforms with "exact" waveforms 
obtained by numerical relativity approaches. The first step of 
the methodology consisted in studying the small-mass-ratio 
limit, v = rnivtii I M 2 <C 1, in which one can use the well 
controllable "laboratory" of numerical simulations of test par- 
ticles (with an added radiation-reaction force) moving in black 
hole backgrounds. Historically, this "laboratory" has been im- 
portant in understanding/discovering several key features of 
GW emission near black holes. A notable example of this 
being the work of Davis, Ruffini and Tiomno [53] which dis- 
covered the transition between the plunge signal and a ringing 
tail when a particle falls into a black hole. The recent study 
of inspiralling and merging small-mass-ratio systems [22] led 
to introducing (and testing) the following improvements in 
EOB dynamics and waveforms: (i) an improved analytical ex- 
pression for the ((£, m) = (2, 2) even-parity Zerilli-Moncrief) 

waveform ^22 which includes a resummation of the tail ef- 
fects, and a 3 +2 PN-accurate "non-linear" amplitude correc- 
tion, (ii) the inclusion of non-quasi-circular corrections to 
the waveform, (iii) the inclusion of non-quasi-circular correc- 
tions to radiation reaction, and (iv) an improved treatment of 
the matching between the plunge and ring-down waveforms 
which takes into account a new understanding of the impor- 
tance of the number of quasi-normal-modes (QNMs), the sign 
of their frequencies, and the length of the interval on which 
the matching is done. The resulting improved implementa- 
tion (when v <C 1) of the EOB approach yielded very faith- 
ful waveforms whose amplitude and phase agreed remarkably 
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TABLE I: Initial ADM mass (scaled by M = mi + 7712) and angular momentum of the spacetime (scaled by M 2 ); final mass (scaled by M) 
and dimensionless spin parameter j{ = J{ /A/ f 2 of the merged black hole; dominant (quasi-normal-mode) complex frequency of the ringdown; 
for two different grid spacings h. 



h/M 


M ADM /M 


Jadm/M 2 


M} OT /M 


•hoi 
Jf 


M f ling /M 


■ ring- 


Mo 2220 


0.024 
0.020 


0.990484 
0.990484 


0.991803 
0.991803 


0.951531 
0.951611 


0.687142 
0.686916 


0.959165 


0.684639 


0.085475 +i 0.551040 



well with the "exact" ones: in particular, the EOB phasing dif- 
fered from the "exact" one by less than ±1.1% of a cycle over 
the whole process. 

The program initiated in [22] was pursued in [24] where the 
comparable-mass version of the improved, resummed 3 +2 -PN 
accurate waveform was compared with the recently published 
inspiral simulation of the Caltech-Cornell group [18]. It was 
found that, by exploiting the combined flexibility in 05 and 
Wpoic, one could reach a remarkable phase agreement, better 
than 0.001 GW cycles over 30 GW cycles. Here, we shall 
similarly exploit the flexibility in 05 and v po \ e to best fit the 
AEI merger waveform. 

Let us recall that the EOB approach is a non-perturbatively 
resummed analytic technique which consists of several differ- 
ent elements: 

• a Hamiltonian H Ica \ describing the conservative part of 
the relative two-body dynamics. The key ingredient of 
this Hamiltonian (defined in Eqs. (13) and (14) of [24]) 
is the "radial potential" A(r). 6 This radial potential is 
defined, at n-Post-Newtonian (PN) order, as the (l,n) 
Pade resummation [4] of its Taylor (i.e. usual PN) ex- 
pansion (written in Eq. (15) of [24]). 

• a radiation reaction force T v (denoted after its 
rescaling by l//i), which is defined as a Pade re- 
summation [1] of its Taylor expansion. See Eq. (17) 
of [24] where /dis is the P4 Pade resummation of 
(1 - w/wpoio)^ Taylor (w; v). The coefficients of F Ts * lox 
in Eq. (18) of [24] have been derived in Refs. [54-59]. 
We shall also consider at the end, following Ref. [22], 
the possibility of modifying T v by a non-quasi-circular 
correcting factor, Eq. (13). 

• improved "post-post-circular" dynamical initial data 
(positions and momenta) as advocated in Sec. Ill B 
of [24]. To explain the improved construction of initial 
data let us introduce a formal book-keeping parameter 
e (to be set to 1 at the end) in front of the radiation reac- 
tion jF v in the EOB equations of motion. One can then 



' Except when said otherwise, we henceforth systematically scale dimen- 
sionful quantities by means of the total rest mass M = mi + m,2 of the 
binary system. For instance, we use the dimensionless EOB radial coor- 
dinate r = Reob/M, with G = 1. Note also that v = fi/M with 
fi = m\m,2 /M. 



show that the quasi-circular inspiralling solution of the 
EOB equations of motion formally satisfies 



Pv =3o(r)+e 2 j 2 (r)+0(e 4 ), 
Pr , =sir 1 (r)+e 3 Tr 3 (r)+0(s 5 ). 



(1) 
(2) 



Here, jo[r) is the usual circular approximation to the 
inspiralling angular momentum as explicitly given by 
Eq. (4.5) of [3], while the order e ("post-circular") term 
7Ti(r) is obtained by: (i) inserting the circular approx- 
imation p v = jo(r) on the left-hand side (l.h.s) of 
Eq. (10) of [21], (ii) using the chain rule djo(r)/dt = 
(djo(r)/dr)(dr/dt), (iii) replacing dr/dt by the right- 
hand side (r.h.s) of Eq. (9) of [21] and (iv) solving for 
p rt at the first order in e. This leads to an explicit result 
of the form (using the notation defined in Ref. [21]) 



ewi(r) 



VHHr. 



1/2 



djo 
dr 



(3) 



where the subscript indicates that the r.h.s. is eval- 
uated at the leading circular approximation e — > 0. 
The post-circular EOB approximation (70, 7Ti) was in- 
troduced in Ref. [3] and then used in most of the subse- 
quent EOB papers [6, 19, 21, 23, 23, 26]. The post- 
post-circular approximation (order e 2 ), introduced in 
Ref. [24] and used here, consists of: (i) formally solv- 
ing Eq. (11) of [21] with respect to the explicit p 2 , 
appearing on the r.h.s., (ii) replacing p rt by its post- 
circular approximation (3), (iii) using the chain rule 
diTi(r)/dt = (diT\(r)/dr) (dr/dt), and (iv) replacing 
dr/dt in terms of 7Ti (to leading order) by using Eq. (9) 
of [21]. The result yields an explicit expression of the 
typep 2 ~ jQ(r)[l + e 2 /c2(»")] of which one finally takes 
the square root. In principle, this procedure can be iter- 
ated to get initial data at any order in e. We found that 
the post-post-circular initial data (joVl + £^2, tti ) are 
sufficient to lead to negligible eccentricity when starting 
the integration of the EOB equations of motion at radius 
r = 15. 

• an improved, resummed "inspiral-plus-plunge" (here- 
after abbreviated as "insplunge") waveform 1 of the 
form 



Here, as before, we work with a metric-level ("h"), rather than 
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{4m) *22 SplUnSC W - -4^(r^) 2 / 2 N 2 QC F 22 e- 2i * , 

(4) 

where $(t) is the EOB orbital phase, f2 = $ is the EOB 
orbital frequency, r w = rijj 1 ^ is a modified EOB ra- 
dius, with ip being defined in Eq. (22) of Ref. [25]. The 
factor F22 is a resummed, 3 +2 -PN-accurate complex 
amplitude correction valid during the (adiabatic) inspi- 
ral, and is an extra complex correcting factor, 

aimed at taking care (in an effective way) of various non 
quasi-circular (NQC) effects during the plunge. F22 
is defined in Eqs. (5)-(ll) of [24], with / 22 being the 
(3,2) Pade resummation of Z^ 7 ' 01 [see also Ref [60] for 
an independent derivation of the nonresummed, 3 PN- 
accurate (2, 2) waveform]. 

• a ringdown waveform 

N N 

where the label N actually refers to a set of in- 
dices (£,£', m, n), with (I, m) = (2,2) being the 
Schwarzschild-background multipolarity degrees of the 
considered (Zerilli-Moncrief-type) waveform ~ 
him, with n = 0,1,2,... being the "overtone num- 
ber" of the considered Kerr-background Quasi-Normal 
Mode (QNM; n — denoting the fundamental mode), 
and !! the degree of its associated spheroidal harmonics 
Se> m (ao~, 0). In addition er^ = ± iu>^ refers to the 
positive/negative complex QNM frequencies (a^ > 
and lu^ > indicate the inverse damping time and the 
oscillation frequency of each mode respectively). The 
sum over £' comes from the fact that an ordinary spheri- 
cal harmonics Yi m (9, </>) (used as expansion basis to de- 
fine $£,„ ) can be expanded in the spheroidal harmonics 
Sii m {aa, 9)e lm ^ characterizing the angular dependence 
of the Kerr-background QNMs [61]. 

• an improved way of matching the inspiral-plus-plunge 
waveform to the ring-down one, on a (2p + l)-tooth 
"comb" (t m - pS, t m - (p - 1)5, . . . , t m - S, t m , t m + 
6, . . . ,t m + pS), of total length At = 2pS, which is 
centered around some "matching" time t m . Below we 
will fix the integer p to the value p = 2, corresponding 
to five matching points. 



• Finally, we define our complete EOB matched wave- 
form (from t = —00 to t = +00) as 

*f 2 OB (i) = 9{t m - i)^ plunge (t) 

+ 6{t - t m )^ do ^(t) (6) 

where 6(t) denotes Heaviside's step function. Note that, 
if one wanted to have a C°° transition between the two 
waveforms one could replace 9(t — t m ) by one of Lau- 
rent Schwartz's well-known smoothed step functions 
(or "partitions of unity") 9 £ ((t — t m )/(2p6)). 

Let us now state the specific choices made here for the var- 
ious EOB ingredients just recalled. Some of these choices 
correspond to various ways of "flexing" the EOB formalism 
(in the sense of Ref. [27]). 

• We "flex" the currently known 3 PN-accurate EOB 
Hamiltonian [4, 45] by introducing an (effective) 4 PN 
Hamiltonian parameter 05, parametrizing an additional 
contribution +a^v/r 5 in the main EOB radial func- 
tion A(r). This parameter has already been intro- 
duced (under varying notations) in several previous 
works [5,23, 24, 27, 28]. 

• Similarly, the EOB radiation reaction force (defined by 
Eq. (17) of [24]) is "flexed" by allowing the Pade- 
resummation parameter v po \ e to differ from the "stan- 
dard" value v^oicM advocated in [1]. 

In addition, we shall also briefly explore another phys- 
ically natural flexibility in the radiation reaction, which 
was introduced (and shown to be physically needed for 
faithfulness) in [22]: the multiplication of the radiation 
reaction by a non quasi-circular (NQC) correction fac- 
tor /^ C , see Eq. (13) below. 

• To define precisely the "insplunge waveform" (4) we 
need to specify: 

(i) the argument x(t) used in the / 22 "brick" within F 22 
(see Eq.(10) of [24]). We shall use here x = ft 2 / 3 
where VL is the EOB orbital frequency . 

(ii) the Pade resummation of the Taylor expansion 
/ 22 ayl ° r °f /22- As in [24] we shall use a P 2 Pade. 

(iii) the definition of the non quasi-circular (NQC) cor- 
rection factor /^ C . To do this we follow the rationale 
explained in [22]. For convenience, we choose (as sug- 
gested in footnote 9 of [22]) afactorized complex NQC 
factor 



curvature-level ("1/14"), waveform. However, we normalize here this 
metric waveform in the same "Zerilli-Moncrief" way as in the test- 
mass work [22]. This differs simply by a numerical factor from 
both the usual tensor-spherical harmonics (I, m) metric amplitude h tm 
and the related metric variables Qt' x extracted from the NR evolu- 



tion [31]: Rhi„ 



1 

V2 



= ^(t + 2)(e + i)i(t - 1) (*w + 



/■NQC 
J22 



Pr, 



(ray 



exp 



(7) 



in which a (denoted a' in the cited footnote) affects only 
the modulus, and b (alias b') only the phase. To ease 
some technical problems during the ring-down linked 
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to the fact that fl(t) tends exponentially towards zero as 
t — * +00 we have added a ("cut-off") constant e to the 
first denominator (rVl)' 2 . As discussed in [22], one can 
a priori analytically determine a "good" value of the 
NQC-modulus parameter a by requiring that the mod- 
ulus of the full EOB insplunge waveform (4) be max- 
imum at the "EOB-light-ring", i.e. when the EOB or- 
bital frequency il reaches a maximum. Ref. [22] men- 
tioned that, in the v -C 1 limit, this requirement im- 
plied a = 1/2 (when e = 0). We found, by numer- 
ically exploring the modulus of <I' 2 I 2 l!plunse (i), that the 
same value, a = 1/2 (together with e = 0.12), can 
be used in the case v = 1/4 considered here. Con- 
cerning the NQC-phase parameter b we simply choose 
b = 0. [Note that the comparable-mass resummed EOB 
waveform of [24] uses a refined estimate for the addi- 
tional phase 622 of & 1 J^ P unBe (t) compared to the one 
used in [22].] 

• Concerning the choice of QNMs we recall that the dis- 
cussion of the physical excitation of QNMs in [22] (see 
the summary in Fig. 4 there) suggested that it is suffi- 
cient to use only positive-frequency QNMs in the ring- 
down waveform (5). This is what we shall do here as 
well. 

A new feature of the comparable-mass case (w.r.t. the 
small v limit) is the "mixing" between various £' QNMs 
(with £' 7^ £) that can enter a given (I, m) multipolar 
wave. This mixing is due to the "aui coupling" terms 
in the separated Teukolsky equations and has been dis- 
cussed in [19, 61]. However, as emphasized in [19], this 
coupling has only a small effect on the (£, m) = (2, 2) 
waveform. We shall neglect it and consider only the 
(positive-frequency) QNM modes having the same val- 
ues of (£, m) as the considered multipolar waveform 
hi m 0-6. (2, 2) in the present paper). 

On the other hand, contrary to other recent implemen- 
tations of the EOB approach [19, 20, 23], we shall 
use a matching comb with^zve teeth (p = 2) and^ve 
(positive-frequency) QNMs a+ mn = aj mn + iw+ m , 
with I = 2, m = +2, and n = 0, 1, 2, 3, 4. To estimate 
the values (as functions of the mass and spin of the final 
black hole) of the damping time and the oscillation fre- 
quency of each mode we did the following: (i) for the 
first three modes we used the approximate fitting formu- 
las given in Appendix E of Ref. [52]; while, (ii) for the 
fourth and fifth modes (i.e. n = 3, 4) we noticed that the 
graphic results of [62] (notably his most relevant Fig. 4) 
exhibit an approximate linearity of the complex QNM 
frequency crf mn as a function of the overtone number 
n. [Indeed, the corresponding points in the complex a 
plane are approximately aligned.] We then exploited 
this approximate linearity to express the needed n = 3 
and n — 4 complex frequencies as linear combinations 
of the above-discussed n = 1 and n = 2 ones. 

• Concerning the matching, on a multi-toothed "comb", 
of the inspiral-plus-plunge waveform to the ring-down 



one we need to specify the two parameters defining such 
a comb, namely the central"matching" time t m , and the 
spacing between the teeth of the comb: S = Af/4 8 . In 
conformity with the basic idea proposed in the original 
EOB paper [3] we choose as central matching time t m 
the so-called "EOB light-ring crossing" time; i.e., the 
EOB dynamical time when the EOB orbital frequency 
ft reaches its maximum. See [22] for a detailed discus- 
sion of why such a choice is physically preferred. Con- 
cerning the choice of the comb spacing S, we expect 
from [22] that a value of order 5 = (7.2M)/4 = 1.8M 
will be good. Below, we shall explore values near this 
one. 



IV. COMPARING THE NR WAVEFORM TO EOB ONES 

As explained in Sec. II, the basic NR data that we shall 
consider is a time-series giving the quadrupolar [(£,m) = 
(2, 2), Zerilli-Moncrief-normalized] metric waveform l f , ^ r 2 R as 
a function of the NR time variable 9 £nr (measured in units of 
M = mi + m,2i- ^^(^nr) is a complex number. The NR 
results consist of the real and imaginary parts of ^\^2~ ■ ^ i s ' 
however, more convenient to decompose the complex wave- 
form in modulus (or amplitude) and phase, say 



*& R (*nr) 



^ 2 R (t N R)cxp(-i0N R (t NR )) 



(8) 



The 2tt ambiguity in the phase is fixed by starting with the 
principal value of the argument of at the beginning of 
the NR simulation, and then keeping track of the 2ir turns as 
the waveform continuously unfolds. 

One can then compute the gravitational wave (GW) fre- 
quency as a function of time by (numerically) differentiating 
the GW phase 



NR/ . 
J22 l^R, 



dt 



(9) 



NR 



[It can equivalently be obtained by computing the imaginary 
part of the logarithmic time derivative of $™(*nr)J 

As emphasized in [24], another useful diagnos- 
tics of GW radiation is the GW phase acceleration 
a = dui/dt = d 2 (f)/dt 2 considered as a function of the 
GW frequency uj. However, because of the presence of 
some additional high-frequency wiggles in <\> and to in the 
NR data, we shall not consider here the phase-acceleration 
curve a(ui). Instead, we shall directly compare the numerical 
GW amplitude, phase and frequency to their analytical, EOB 
counterparts. 



8 Note that in [22] we used the letter <5 to denote the full width At of the 
comb. 

9 As mentioned in Sec. II, we use the waveform extracted at a (coordinate) 
radius -Rnr = 120Af c ~ 120M, and tNR is the time of the "observer" 
located at the latter radius. 
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The integration of the basic EOB dynamical equations 
(written in [24]) gives, for each chosen value of the EOB 
"flexibility parameters" (notably a.5 and u po i e )> several impor- 
tant time series, and notably: (i) the EOB orbital frequency 
S1(£eob), where £eob is the EOB dynamical time scale (mea- 
sured in units of M); (ii) the new, resummed matched 3 +2 - 
PN-accurate quadrupolar EOB waveform ^f^^EOB); then, 
from the latter, one can define (as for the NR case) the 
corresponding EOB amplitude, Af 2 OB (iEOB), EOB phase, 
</ , f 2 ° B ( < EOB), and EOB frequency w e 2 ob (>eob)- To com- 
pare the NR and EOB phase time-series 022 R (*nr) an d 
( / , 22° B (^eob) one needs to shift, by additive constants, both 
one of the time variables, and one of the phases. In other 
words, we need to determine r and a such that the "shifted" 
EOB quantities 

*bob = *eob + t , <t>'i OB = 0! OB + a (10) 

"best fit" the NR ones. One convenient way to do so is first 
to "pinch" the EOB/NR phase difference at two different in- 
stants (corresponding to two different frequencies). More pre- 
cisely, one can choose two NR times t± R , t^ K , which deter- 
mine two corresponding GW frequencies 10 u)\ = W22 (*i )> 
u>2 = cJ22 a (#2 )> an d then find the time shift t(u>i,U2) 
such that the shifted EOB phase difference, between uj\ 

and wa, A0 EOB (r) = ^ OB (4 EOB ) - 2 ! OB (i'i EOB ) = 
^ E 2 OB (t EOB +r)-(/) E 2 OB (t EOB +r) is equal to the correspond- 
ing (unshifted) NR phase difference A(/) NR = </^ 2 R (^ R ) - 
022 (if R ) . This yields one equation for one unknown (r), 
and (uniquely) determines a value t(u)±, 0J2) of r. [Note that 
the u>2 — > wi = u! m limit of this procedure yields the one- 
frequency matching procedure used in [18].] After having so 
determined t, one can uniquely define a corresponding best-fit 
phase shift a(wi, w 2 ) by requiring that, say, 22 OB (ii EOB ) = 

E 2 ° B (4 EOB )+« = ^2 R m- 

Having so related the EOB time and phase variables to the 
NR ones we can straigthforwardly compare all the EOB time 
series to their NR correspondants. In particular, we can com- 
pute the (shifted) EOB-NR phase difference 

A-.-4 OBNa M = ^ EOB (t' EOB ) - <^ R (; NR ). (ii) 

In the following we will chose two matching instants (and 
corresponding frequencies) that take place during late inspiral 
and plunge, namely: t^ R = 999.72, t^ R = 1494.94 corre- 
sponding to uj\ = 0.06815, ll>2 = 0.2457 (all expressed in M 
units). 

To numerically implement the EOB/NR comparison we 
need to choose some values for the various "flexibility param- 
eters" of the EOB framework. We have summarized above 
what are these parameters, and we have already indicated the 
values we chose for some of them. Among the remaining ones 
that need to be chosen, the two most crucial ones are 05 and 



Wpoio- Recently, Damour and Nagar have shown, by using 
some of the data published in [18], that the inspiral waveform 
(for GW frequencies smaller than about 0.14/M) could be re- 
markably well matched by the EOB one if one chose values 
of 05 and Wp ic following the rather precise correlation plot- 
ted in the upper panel of Fig. 3 in [24]. Here, as we are ex- 
ploring a different physical regime (late inspiral, plunge and 
coalescence, with GW frequencies mostly larger than about 
0.1/M), and comparing to a different set of numerical data, 
we shall not a priori impose the precise correlation between 05 
and Wpoio found in [24]. However, we shall make use of some 
previous results suggesting a preferred range for the values 
of a.5. On the one hand, Ref. [23] showed that the faithful- 
ness (in the sense of Sec. VIA of [27]) of (restricted) EOB 
waveforms against NASA-Goddard NR coalescence wave- 
forms was largest when 05 belonged to some rather wide inter- 
val (which also depended on the considered mass ratio). See 
Fig. 2 (right panel) in [23] from which one might conclude 
that 05 lies probably between ~ 10 and ~ 100. Buonanno et 
al. then chose 05 = 60 as "best fit" value. On the other hand, 
Ref. [24] found that the phase agreement between (resummed) 
EOB waveforms and a rather long inspiral NR waveform was 
at its best when 05 lied in a similarly wide interval (between 
~ 10 and ~ 80) centered around 05 ~ 40. In view of these 
results we shall focus, in the following, on two representative 
values of a.5, namely 05 = 25 (representative of the leftward- 
side of preferred 05 values), and a.5 = 60 (representative of 
the rightward-side of preferred 05 values, and chosen as best 
value by [23]). We have also checked that the values of 05 
between 25 and 60 lead (with appropriate choice of u po i c ) to 
results that are at least as good as the ones we shall exhibit 
below. 



A. Comparing NR to resummed EOB for a 5 = 25 

At this stage we have essentially fixed all the flexibility 
of the EOB formalism apart from the choices of u po i e , and 
of the comb spacing 5. Among these two parameters, only 
the former one, v po i e , is important for getting a very accurate 
phase agreement between EOB and NR. When 05 = 25, we 
found (by trial and error) that 11 u po ie = 0.6241 (together with 
S = 1.7 AI{ which is, however, less crucial) yields an excellent 
EOB/NR agreement. We exhibit our results in the four panels 
of Fig. 1. 

The top-left panel of Fig. 1 compares the NR GW fre- 
quency both to the (matched) EOB GW frequency, and to 
twice the orbital frequency. The time axis is t^R, and/or 
(see above) t^oB = *eob + t (with r = — 2032M for the 
present case). The vertical lines on the right indicate the cen- 
ter and the outlying "teeth" of our matching comb, which is, 



Alternatively, one can start by giving oneself u>\ , u>2 and determine the NR 
instants t ^ R , t ^ R at which they are reached. 



Though we did not investigate thoroughly what "error bar" can be put on 
such a "best" value of ^ po i e , the numerical studies we did indicate that a 
change of ±2 on the last (i.e. fourth) digit that we quote is sufficient to 
entail a visible worsening of the phase difference A02^- IBNH . 
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FIG. 1: Comparison between EOB and NR waveforms for a$ = 25 and v po \ c = 0.6241: frequencies (top-left), phase difference (top-right), 
amplitudes (bottom-left) and real parts (bottom-right) of the two gravitational waveforms. The vertical line at £nr = 1509 locates the 
maximum of (twice) the orbital frequency Q. (alias the "EOB-light-ring") and indicates the center of our matching comb (whose total width is 
indicated by the two neighboring vertical lines in the top-left panel). The vertical dashed line at t NR = 1482 indicates the crossing time of the 
adiabatic LSO orbital frequency (f^LSO = 0.1003). 



as explained above, centered on the maximum of the EOB or- 
bital frequency (also called "EOB-light-ring"). The interval 
between the two vertical lines (LSO and "EOB-light-ring") 
defines the "plunge". The dashed vertical line on the left (at 
^nr = 1482) indicates the crossing time of the adiabatic Last- 
Stable-Orbit (w-LSO in the sense of [3]). Note that the three 
frequencies are initially close to each other, but that, later, 
2il separates from w^ 2 R and wf 2 ° B ' which continue to be in 
very close agreement, except for a slight discrepancy around 
merger, which, within the EOB approach, is conventionally 
supposed to take place at the maximum of H. Note also the 
good agreement between the EOB GW frequency during the 
ringdown plateau, and the average of the NR one. As dis- 
cussed in Sec. II the values for the mass and dimensionless 
spin of the final black hole that we used (together with [52]) 
to compute the QNMs frequencies are: M f ring = 0.959165M, 



j f rin s = 0.684639. 

The top-right panel of Fig. 1 shows the EOB-NR phase dif- 
ference, Eq. (11), ("pinched" at the two instants, tf R , t^ R , 
given above). It is remarkable that the (two-sided 12 ) EOB- 
NR phase difference over the time interval (639M, 1524M) 
(which covers about 12 GW cycles of inspiral, plunge, and 
early ring-down) is smaller than about ±^0.068 radians, 
which corresponds to ±0.005 GW cycles. 

The bottom-left panel of Fig. 1 compares the NR GW 
amplitude to the resummed 3 +2 -PN accurate EOB one. It 



As the reference level of any phase difference A.cj> is arbitrary, it is conve- 
nient to use a "middle" reference level such that A<^(i) varies between — e 
and +e over the considered interval. Werefertoie = ±l/2[max(A</<) — 
min(A</>)] as the two-sided phase difference. 
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also shows the orbital frequency fl as an aid to locate the 
merger. One notices a very good agreement between the two 
amplitudes. During the interval (1100M, 1400M) the frac- 
tional EOB-NR amplitude difference varies between — 1% 
and +1%. After i NR = 1400M, this fractional difference 
increases from +1% to a maximum of +18% (reached at 
iNR — 1509M) and then decreases to take values of order 
—5% during the observationally relevant part of the ringdown. 
Note also that the NR equal-mass amplitude (divided by v, i.e. 
by \x) time series is qualitatively, and even quantitatively, very 
similar to the corresponding NR test-mass amplitude time se- 
ries shown in Fig. 3 of [22] . For instance, the value of the max- 
imum amplitude is ~ 0.3 in both cases. A similar qualitative, 
but not quantitative, parallelism exists for the two correspond- 
ing frequency time series (the v = 1/4 frequency levelling off 
at a higher "plateau"). 

Finally, the bottom-right panel of Fig. 1 compares the real 
parts of the NR and EOB waveforms. The two vertical 
lines delimit the interval between LSO and "EOB-light-ring". 
Again the agreement between the two waveforms is impres- 
sive. Note that this last panel shows only the late inspiral, 
plunge and ringdown. From the panel showing the phase dif- 
ference, one can gather that the agreement stays as impressive 
over a much longer time span of order 1000M (essentially 
from iNR ~ 500M to the end of ringdown). 



B. Comparing NR to resummed EOB for a 5 = 60 

Let us now consider our second representative value of the 
effective 4 PN radial potential parameter, as = 60. As be- 
fore we chose 6 = 1.7M{, We also selected the same phase 
"pinching" interval as above. Then, by trial and error, we 
found that iJ po ic = 0.5356 yields an excellent EOB/NR agree- 
ment 13 . 

We exhibit our results in the four panels of Fig. 2, which 
are entirely parallel to those of Fig. 1 . The remarkable level of 
EOB/NR agreement that we get now, when a$ = 60, is rather 
close to the one that we got above when as = 25. At this 
stage, there is no rationale for saying that either value of 05 
is preferred over the other (though as = 25 yields somewhat 
better results). Some partial numerical tests that we performed 
suggest that this conclusion extends to (at least) all values of 
as between 25 and 60. 

Some of the numbers quantifying the EOB/NR agreement 
are: 

(i) the (two-sided) EOB-NR phase difference over the time 
interval (500M, 1550M) (which covers about 13 GW cycles 
of inspiral, plunge, and most of the ring-down) is smaller than 
about ±|0.13 radians, which corresponds to ±0.01 GW cy- 
cles; 

(ii) during the interval (1100M, 1400M) the fractional 



Note that this "best" value of v po \ e (f° r a 5 = 60 and u = 1/4) happens 
to be numerically close to the best fitting v po i c ~ 0.53 value that Ref. [24] 
found in the test-mass limit v — ► 0. 



EOB-NR amplitude difference varies between —0.8% and 
+0.55%. After f NR = 1400M, this fractional difference in- 
creases from +0.55% to a maximum of +23% (reached at 
iNR — 151 IM) and then decreases to take values of order 
+6% during the observationally relevant part of the ringdown. 



C. Contrasting resummed EOB with restricted EOB, for 
a 5 = 60, by comparing NR to a standard restricted EOB 
waveform 

Finally, we wish to illustrate the importance (for reaching 
a high level of accuracy) of the various ingredients used in 
our present, resummed version of EOB (using a time-extended 
"comb matching" to 5 QNMs) by comparing NR to the type of 
simpler implementation of the EOB framework used in [23]. 
Using again as = 60 (which was chosen as best value in [23]), 
we compare NR to the following implementation of EOB: 

• we use for v po \ e the "standard" value v™\ c (i>) advo- 
cated in [1]. 

• we use (as originally proposed in Ref. [3]) the following 
(Newtonian-order and Kepler-law-assuming) restricted 
quadrupole waveform 

^ 2 K (t) = -AvJ^f' 3 cxp(-2i$) , (12) 

without any explicit PN (F22) corrections, nor any NQC 
(a, b) corrections. 

• we use only 3 (positive-frequency) QNMs. 

• and, we match the plunge and ring-down waveforms 
in a very small interval {8 /Mi = 0.2 instead of our 
preferred 1.7) around the maximum of the orbital fre- 
quency. [Indeed, the matching of the two waveforms 
and their derivatives at a sharply defined moment is 
equivalent to considering the S — > limit of our comb- 
matching technique]. 

The results of such a coarser EOB implementation are 
shown in Fig. 3 (which is parallel to the previous two figures). 
By contrasting Fig. 3 with Fig. 2 (which used the same value 
of as), we see that: 

• the EOB frequency agrees less well with the NR one 
than before, especially around the matching point. 
Note in particular that the post-matching analytical fre- 
quency jumps up from the maximum (doubled) or- 
bital frequency significantly more vertically than be- 
fore, thereby decoupling too soon from the exact fre- 
quency, and accruing a larger dephasing than before 
(because of the too localized matching, and - to a lesser 
degree - the use of only 3 QNMs). 

• the EOB-NR (maximal) phase difference over the same 
time interval (500M, 1524M) is about 2.2 times larger 
than before. One now ends up with a phase difference 
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FIG. 2: Comparison between EOB and NR waveforms for as = 60, u po ic = 0.5356: frequencies (top-left), phase difference (top-right), 
amplitudes (bottom-left) and real parts (bottom-right) of the two gravitational waveforms. The vertical line at £nr = 1510 locates the 
maximum of (twice) the orbital frequency (alias the "EOB-light-ring") and indicates the center of our matching comb (whose total width is 
indicated by the two neighboring vertical lines in the top-left panel). The vertical dashed line at t NR = 1487 indicates the crossing time of the 
adiabatic LSO orbital frequency (f^LSO = 0.1081). 



of ±^0.29 radians, i.e. 0.023 GW cycles over about 
13 GW cycles. The top-right panel of Fig. 3 illustrates 
the fact that matching with 5 QNMs (dashed line) re- 
duces the dephasing accumulated during the transition 
from merger to ringdown. 

• the modulus of the analytical waveform is now dis- 
tinctly larger than the NR one during the inspiral (be- 
cause of the lack of PN corrections). 

• the modulus also exhibits a more significant discrep- 
ancy (+35%) with the NR one at the end of the plunge 
(because of the use of the Kepler-law-assuming oc 
51 2 / 3 , which, as pointed out in [25], tends to overesti- 
mate the amplitude). 

• Note also that one visually notices these differences at 



the level of the GW waveforms. 

• The same resummed-EOB/restricted-EOB comparison 
was done in [22], in the v <C 1 case, with similar con- 
clusions. 

In spite of these relative blemishes, note, however that this 
"coarser" EOB-type implementation still succeeds in follow- 
ing the phase of the exact signal to ±0.023 GW cycles over 
about 13 GW cycles. 

Note that the corresponding EOB/NR agreement exhibited 
in Fig. 4 of Ref. [23] seems to be somewhat better 14 than the 



The reader should however keep in mind that in Fig. 4 of Ref. [23] the 
EOB-NR phase difference is divided by 27r compared to the one showed in 
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FIG. 3: Comparison between the EOB restricted waveform approximation, Eq. (12), and NR for as = 60 and v po i c = v^^ e {v — 1/4) = 
0.6907: frequencies (top-left), phase difference (top-right), amplitudes (bottom-left) and real parts (bottom-right) of the two gravitational 
waveforms. The vertical line at 4nr = 1510 locates the maximum of (twice) the orbital frequency SI (alias the "EOB-light-ring") and 
indicates the matching time. The vertical dashed line at f NR = 1490 indicates the crossing time of the adiabatic LSO orbital frequency 
($7 LS o = 0.1081). 



one exhibited by our Fig. 3. This difference might have sev- 
eral origins, notably: (i) a difference in the accuracy of the NR 
data 15 , and (ii) a difference in the procedure used to best shift 
time and phase between EOB and NR data. 



our Fig. 3. 

15 The data used in [23] did not benefit from the reduction in eccentricity used 
in the data considered here. 



V. CONCLUSIONS 

We have compared a recently proposed, resummed 3 +2 -PN 
accurate Effective-One-Body (EOB) waveform to the result 
of a numerical simulation of a coalescing equal-mass binary 
black hole performed at the Albert Einstein Institute. We find 
a remarkable agreement, both in phase and in amplitude, be- 
tween the new EOB waveform and the numerical data. More 
precisely, we find that the maximal dephasing between EOB 
and numerical relativity (NR) can be reduced below ±0.005 
GW cycles over the last ~ 900M (corresponding to about 
12 GW cycles plus ringdown ones) of the simulation. This 
level of agreement was exhibited for two representative val- 
ues of the effective 4 PN parameter 05, namely 05 = 25 
and 05 = 60, and for a corresponding, appropriately "flexed" 
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value of the radiation-reaction resummation parameter v po i c . 
In addition, our resummed EOB amplitude agrees to better 
than the 1% level with the NR one up to the late inspiral. 

We have also compared the NR data to a coarser implemen- 
tation of the EOB approach (restricted waveform, standard 
Vp*i e , instantaneous matching to 3 QNMs). The EOB/NR 
agreement is slightly less good in this case, though the phase 
agreement remains quite good (±0.023 GW cycles over the 
last — 1000M of the simulation). 

Let us point out a notable feature of our results. In 
the recent work of Damour and Nagar [24], the same re- 
summed 3 +2 -PN accurate EOB waveform was compared to 
a long, very accurate equal-mass inspiral simulation of the 
Cal tech-Cornell group [18]. It was found that an excellent 
EOB/NR agreement was obtained when 05 and u po i e were 
following the rather precise correlation plotted in the up- 
per panel of Fig. 3 of Ref. [24]. Let us denote this cor- 



relation as 0,5 



v 



best inspiral 
pole 



(05). In the present paper, 



we similarly found that the EOB/NR agreement was at its 
best when, for a given as, 16 u po i c was taking a rather pre- 

i- tti_ ^ c± 1 " best insplunge/ \ 

cise corresponding best fit value , say v le ( a s)- 
In particular, we found b-tins P iunge (25) 



best insplungc 
pole 



(60) 



po 

0.5356. 



results of [24] yield inspiral (25) 



best inspiral 
pole 



(60) 



sets of values are v pole 
0.0901 and 



0.6241 and 
On the other hand, the 
0.5340, and 
0.4856. The differences between these 

best insplungc / or \ best inspiral / 



(25) 



pole 

(60) - t;^ e tinspiral (60) = 0.0500. 
Note also that the "best insplunge" v po i c values are in between 
the "best inspiral" ones and the originally advocated [1] one 
v™f e (v = 1/4) = 0.6907. This finding will deserve fur- 
ther investigation in the future. At this stage we can only 
speculate on the various possible origins of this difference: 
(i) it might be due to the fact that, not having access to the 
original NR data of [18], Damour and Nagar had to rely 
on rather coarse measurements extracted from published fig- 
ures; (ii) it might be due to systematic errors in the NR data 
of [18]; (iii) it might alternatively come from systematic er- 
rors in the NR data used in the present paper; (iv) it might 
come from the fact that the "best-fit" ^(upoie) is not a uni- 
form approximation (as a function of frequency) to the ex- 
act radiation reaction (see, in the v — ► limit, the bottom 
panels of Fig. 1 in [24]) and, finally, (v) it might come from 
some "missing physics" in the resummed EOB waveform ex- 
plored here. There are several candidates for this missing 
physics. One suggestion (which follows the original sugges- 
tion of [5]) is that one might need to consider still higher 
(uncalculated) PN contributions to the radial EOB potential 17 



'(25) 



A(u) = 1 — 2u + 2vv? + a^vifi + a^vu 5 + agi^u 6 + • ■ ■ 
where u = 1/r. Another suggestion is that non-quasi-circular 
(NQC) corrections to radiation reaction might modify the 
phasing during late inspiral and plunge. As an example, we 
have looked at this possibility. More precisely, following [22], 
we can introduce a new flexibility parameter a RR 18 such that 
the radiation reaction force is multiplied by a correction factor 

/rr c s iven b y 



/1 



NQC 
RR 



= 1 



-RR 



(rft) 2 



(13) 



Such a factor will be very close to one during the inspiral (and 
therefore will be negligible in the EOB comparison to the 
Caltech-Cornell data), but will start being significantly less 
than one (if a RR > 0) during the late inspiral and plunge, 
which are of interest for the comparison to the presently con- 
sidered data. And indeed, we have found that by choosing 
a value a RR ~ +40 (and e = 0.12 as in the waveform 
NQC factor considered above) we could, when 05 = 25, 
obtain an excellent EOB/NR fit by using the "best inspi- 
ral" value Wpo^ mi! P iral (25) = 0.5340 (instead of the above 

tip°j* msplunsc (25) = 0.6241). This issue needs to be further 
investigated by using the most accurate possible data covering 
both inspiral and plunge. We hope to come back to it in the 
future. 

Finally, we think that the present work, taken in conjunc- 
tion with other recent works on the EOB/NR comparison, con- 
firms the ability of the EOB formalism to accurately capture 
the general-relativistic waveforms. The present work has also 
shown that the recently proposed resummed 3 +2 -PN accu- 
rate waveform is important for defining analytical EOB wave- 
forms that faithfully represent (both in phase and in ampli- 
tude) the waveforms emitted by equal-mass coalescing (non- 
spinning) black hole binaries. 
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Though we did not explore all possible values of 05 , we sampled interme- 
diate values between the representative a§ values we picked and convinced 
ourselves that the same conclusion held for them. 

For simplicity, we consider only linear-in-^ higher PN contributions. If 
the need arises (and the fact that the unequal-mass EOB/NR comparisons 
of [23] seem to exhibit a strong dependence on the mass ratio might suggest 
it) one can easily add in a non-linear v dependence. 



Actually [22] introduced a parameter a RR which is, roughly, the neg- 
ative of a RR , with a NQC radiation reaction factor of the form 1 + 

a nn p 2 r J{rtt) 2 
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